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Abstract 

In a class of beyond-standard-model theories, CP-odd observables, such as the neutron electric 
dipole moment, receive significant contributions from flavor-neutral P-odd and CP-odd four-quark 
operators. However, considerable uncertainties exist in the hadronic matrix elements of these op- 
erators strongly affecting the experimental constraints on CP-violating parameters in the theories. 
Here we study their hadronic matrix elements in combined chiral perturbation theory and nucleon 
models. We first classify the operators in chiral representations and present the leading-order QCD 
evolutions. We then match the four-quark operators to the corresponding ones in chiral hadronic 
theory, finding symmetry relations among the matrix elements. Although this makes lattice QCD 
calculations feasible, we choose to estimate the non-perturbative matching coefficients in simple 
quark models. We finally compare the results for the neutron electric dipole moment and P-odd 
and CP-odd pion-nucleon couplings with the previous studies using naive factorization and QCD 
sum rules. Our study shall provide valuable insights on the present hadronic physics uncertainties 
in these observables. 
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I. INTRODUCTION 



Neutron electric dipole moment (nEDM) has attracted considerable attention over more 
than half a century. For an elementary particle to have non-vanishing intrinsic EDM, simple 
analysis shows that parity-violating as well as time-reversal-violating interactions must be 
present. [T- violation is equivalent to CP- violation (combined charge-conjugation and parity) 
in local quantum field theory] However, in the standard model (SM) of particle physics, such 
interactions arise only from flavor- changing Cabbibo-Kobayashi-Moskawa (CKM) matrix 
elements, which are strongly suppressed phenomenologically, yielding a very small neutron 
EDM of order 10~ 31 ecm. Therefore, an experimental observation of a large-size neutron 
EDM is an unambiguous signal for new physics, widely expected to exit somewhere between 
the electroweak symmetry breaking and TeV scales. 

There has been much speculation in the literature on the nature of new physics. The lead- 
ing theory is supersymmetry (SUSY) introduced to stabilize the electroweak scale against 
the fundamental Planck scale (for a good review, see Ref. [lj]). Various extra dimensional 
models have been popular in recent years as well (for good reviews see Ref. |2j). There 
are other less dramatic possibilities including technicolor jjjj], left-right symmetry [H, min- 
imal extensions of the SM j^], little Higgs models [f|, etc. In each of these models, there 
is new CP violation physics giving rise to a significant neutron EDM. To test CP-violation 
mechanisms, it is important to get accurate predictions of CP-violating observables in these 
models. In particular, one has to deal with the non-perturbative quantum chromodynamics 
(QCD) physics present in the structure of the neutron. 

An efficient way to calculate the neutron EDM is to use the methodology of effective 
field theories (EFT). In this approach, one generates P-odd and CP-odd quark and gluon 
operators after integrating out the heavy particles (including heavy quarks, gauge bosons 
and new particles) and run these operators to a scale around 1 GeV where non-perturbative 
QCD physics becomes important. The effective degrees of freedom involves the light quarks 
(up, down and strange) and gluons. The CP-odd part of the lagrangian is generally written 
as a sum of CP-odd operators of different mechanical dimensions, 

oo 

£ C P-odd = ? J / J C! di{lj)Odi(fJ') > (1) 
d=3 i 

where d = 3,4, 5, etc, is the mechanical dimension of the operators, fi is the renormalization 
scale (taken as ^tcF^ in this paper) and i sums over operators of the same dimension. The 
dim-3 operator is the usual CP-odd quark mass term qij^q, which can be rotated away 
through chiral rotations apart from the £7a(1) anomaly. The dim-4 operator is the usual 9 
term GG. Dim-5 operators include quark electric and chromoelectric dipole operators. Dim- 
6 operators contain various four-quark operators and Weinberg three-gluon operator. The 
matrix elements of dim-4 and, to less extent, dim-5 operators have been studied extensively 
in the literature [7-16], and the uncertainty of the estimates are typically at the level of 



factor 2. The contributions of these operators have also been studied extensively in the 



context of various new physics models (see Refs. [171, [18( for good reviews). 

However, the matrix elements of dim-6 operators have been a challenge to estimate. In 
some beyond-SM theories such as the left-right symmetric model, dim-6 four-quark operators 
dominate the contributions to nEDM. In the literature, the only serious approach that has 
been proposed to calculate their matrix elements is the naive factorization method: breaking 
the four-quark matrix elements into the product of two-quark matrix elements between the 
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nucleoli states and between pion and vacuum While the factorization involving 



mesons can be and has been tested using lattice QCD |22| and the results may be trustable 
to within a factor of 2, the same is not known for matrix elements involving the nucleon 
states. The goal of this paper is to develop a chiral perturbation method combined with 
simple quark models to estimate the four-quark contribution to the nEDM with hopefully 
an improved accuracy. 

The approach we are going to take is the standard chiral perturbation theory (%PT) 



(see, for example, Ref. |23l . |24| ) which has been used to calculate the contribution of 6- 
term to nEDM [13]. One of the successes of the chiral approach can be illustrated by the 
polarizabilities of the nucleon. The electric polarizabilities of the proton and neutron have 



been extracted from experimental data, a e * p = (10.4 ±0.6) x 10 fm , a^ xp = (12.3 ± 1.3) x 
10 _4 fm 3 . The leading contribution in %PT comes from the pion- nucleon intermediate states, 

which diverges linearly as — > and agrees well with the experimental data. One would 
expect then a similar pion dominance in the neutron EDM because the latter also involves 
the intermediate electric dipole excitations. Indeed a pioneering calculation by Crewther 
et al. found that the dominant contribution from the charged-pion chiral-loop diverges 
logarithmically as m n goes to zero, and is proportional to the CP-odd pion-nucleon-nucleon 
coupling g-xNN uM- m this paper, we take this contribution as dominating and consider the 
four-quark operator contribution to g n NN- Of course, there are chiral-regular contributions to 
the nEDM which are of the same order in chiral power counting and numerically competitive 



or even dominating in the real world [18(. We will consider these contributions as well, 
although the mo del- dependence becomes unavoidable. 

In the chiral approach, one first writes down the CP-odd and even lagrangian in terms 
of meson and nucleon fields, 

£ -^Goldstone— boson CP— odd term ^nucleon CP— odd mass term 

+£edm term + ^CP-odd tt-n coupling + (CP — even terms) (3) 

where the Goldstone boson CP-odd lagrangian will generate terms annihilating 7r° and 77 
in the vacuum, or in other words, will produce meson condensates. The condensates will 
turn some of the CP-even terms (as we shall see, those proportional to quark masses) in 
the chiral lagrangian into CP-odd contributions. This will generate an additional CP-odd 
nucleon-mass term, neutron EDM term and CP-odd pion-nucleon coupling. Once this is 
done, one can rotate away the CP-odd nucleon mass term, generating further contributions 
to the neutron EDM terms and the CP-odd pion-nucleon coupling. 

The presentation of the paper is organized as follows: In Sec. II, we classify all flavor- 
neutral P-odd and CP-odd four-quark operators in chiral representations. We also present 
the leading-order QCD scale evolution of these operators. In Sec. Ill, we match these 
operators to the corresponding Goldstone boson operators, baryon operators, and EDM 
operators in ;^PT. We also discuss in the case of Peccei-Quinn symmetry the size of the 
induced 9 term in the presence of these four-quark operators. In Sec. IV, we calculate their 
contributions to the P-odd and CP-odd nucleon-pion vertices and the CP-odd nucleon mass 
using factorization in the case of meson matrix elements and simple quark models for the 
nucleon ones. In Sec. V, we study the four-quark contribution to the neutron EDM in 
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the chiral approach supplemented with factorization and quark model estimates of counter 
terms, and the results are compared with other calculations in the literature. The comparison 
and analysis show that the hadronic physics uncertainties here can be quantified to within 
a factor of two for operators generating unsuppressed meson condensate contributions. We 
conclude the paper in Sec. VI. 

To make the convention clear, the EDM operator for spin-1/2 particles is defined as 

C = - l -d E i,a, u i 1 ^F^ ) (4) 

where = ~ [7 M ,7i/] and 75 = i7°7 1 7 2 7 3 in the standard Dirac representation. 



II. P-ODD AND CP-ODD FOUR-QUARK OPERATORS: CLASSIFICATION, 
RUNNING AND MIXING 

We consider three light quark flavors: up, down and stran ge. Flavor-neutral P-odd and 
CP-odd four-quark operators can be divided into two groups [l9|: The first group includes 
18 operators made of two different flavors: 

On = (qii 5 q)(qq') , 

O12 = (qq)(qil5q) , 

021 = &l 5 t a q)(q't a q') , 

022 = (qt a q)(q'tl 5 t a q') , 

4 = (qt^t a q)(q'a fl Xq') , (5) 

where q, q' = u,d,s and q < q'. The second group includes 6 operators made of a single 
quark flavor 

0[ = (qil5q)(qq) , 

0' 2 = (qil,t a q)(qt a q) , (6) 

where t a are generators of the SU(3) color group. All other flavor-neutral P-odd and CP-odd 
four-quark operators can be related to the above through Fierz transformation [191 ]. 

To match the above quark operators into the hadronic ones in x?T, we have to classify 
the former into irreducible representations of the chiral group SU(3)l x SU(3)r. Take the 
operator ui^udd as an example, which can be decomposed as 

ui^udd = —iuRULdRdi + iuLURdftdL + h.c. , (7) 

where qi,R = Pi,nq with Pl,r = (1 =F 7s)/2. The first term can be further decomposed as 

— % — — 

-iu R u L d R d L = --{u R u L d R d L + d R u L URd L ) 

% — — 

--(u R u L d R d L - d R u L u R d L ) 
= -- A S%qmq k L q R3 q l L ~ U^ n A)qR k qTq R iql , (8) 
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where effi, = e> kl € w 



and 
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with other elements vanishing. The second term of Eq. ([7]) can be written as 



iu L u R d R d L = iH J u H 2k q Lj q R q Rl q l L 



where 
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(10) 
(11) 
(12) 



In this way the operator uijsudd is decomposed into (3, 3), (6, 6) and (8, 8) representations 
of SU (3)^ x SU (3)^, and A, S, H 1 , and H 2 can be regarded as spurion fields in the sense that 
if they transform as (3,3), (6,6) and (8,8) under chiral transformation, the corresponding 
terms in Eqs. ([H]) and fill I) become invariant. These spurion fields will be used in xPT 
to construct the effective operators corresponding to the same four-quark operators. All 
spurion fields for four-quark operators with different Dirac and color structures are shown in 
Table [H [It is easy to see that there is no (1, 1) operator because any such operator must be 
expressible in terms of products of chiral-even quark currents, which cannot yield CP-odd 
contributions.] 

The four-quark operators usually emerge at a high energy scale where some heavy par- 
ticles have been integrated out. To match them to hadronic operators in effective theories, 
one must run them down to a low energy scale where non-perturbative physics becomes 
important. We can choose this to be 1 GeV or the lattice cut-off 1/a, where a is the lattice 
spacing. In this work, we take [i = ^7rF w , with F n = 93 MeV. These operators mix with each 
other when the energy scale changes. Although many of the mixings have been calculated 
in the literature before [see Ref. 25|, for example], we recalculate them and present the 



complete result here for easy reference: 
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Clearly operators with different quark flavor structures do not mix. Since SU(3)l x SU(3) r 
symmetry is broken only by quark masses, four-quark operators belonging to different chiral 
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SE = 4 
other components are zero 



#i = 0, # 2 = 



TABLE I: Spurions for CP-odd 4-quark operators. The first six together with three tensor struc- 
tures yield 18 operators in Eq. (4) and the last three with two tensor structures yield six operators 
in Eq. (5). 



irreducible representations do not mix either. Therefore, we can further simplify Eq. (12), 



/ (3,6) X 





(15) 



(16) 



where Op' 6,8 ^ means the projections of the operator Oi on the representations (3,3), (6,6) 
and (8,8), respectively. It is easy to see that the (3,3) and (6,6) projections of O a and 
Oi2 are the same with i = 1,2, whereas their (8, 8) projections differ only by the sign. The 
tensor operators do not have (8, 8) components and therefore do not participate in Eq. ( jl"6l) . 
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The four-quark operators may also mix with P-odd and CP-odd operators with dimension 
less or equal to 6. For mixing with lower- dimensional operators, either quark masses or power 
divergences will appear. The only other dimension-6 operator is the Weinberg operator [26[ 

1 

6" 



O w = --f abc e^G1 p G b J>G c a p, (17) 



which is a singlet under chiral transformation. Since the four-quark operators contain no sin- 
glet component, the mixing between them and Ow vanishes. The evolution of the Weinberg 



operator can be found in Ref. [27 



The P-odd and CP-odd dimension-5 operators are the quark electric dipole moment 
operators (QEDM) and quark chromo-electric dipole moment operators (QCDM). In prin- 
ciple, they belong to (3, 3) of the chiral group. However, they can mix logarithmically with 



four-quark operators multiplied by the quark mass which transforms also like (3, 3) [28|, |29 

Finally, the four-quark operators can have mixing with mqi'j^q with quadratically diver- 
gent coefficients. Usually, one defines the four-quark operators with quadratic divergences 
subtracted, as is natural in dimensional regularization where all quadratically divergent in- 
tegrals vanish by definition. Equivalently, this can be achieved, for example, by demanding 
the CP-odd four-quark operators have vanishing contribution between QCD vacuum and 
CP-odd meson states in perturbation theory. However, as we shall see in the following sec- 
tion, they can have non-perturbative contributions. The exact physical implication of this 
non-perturbative contribution will be discussed in Sec. V.E. 



III. MATCHING TO OPERATORS IN CHIRAL PERTURBATION THEORY 

Generically, any P-odd, CP-odd quark-gluon operator contributes to all P-odd, CP-odd 
hadronic operators in ;\PT; the latter are constructed in terms of Goldstone-boson (pion, 
kaon, eta) fields and baryon fields. Here we consider just the contributions to the Goldstone- 
boson CP-odd interactions, nucleon CP-odd mass term, n-N CP-odd coupling, as well as 
the neutron EDM term, 



^-'Goldstone-boson CP— odd term -^nucleon CP— odd mass term 
+£cp —odd 7T— N coupling term • 



(18) 



Following the standard practice in the literature, we imbed the Goldstone-boson fields in 
the unitary matrix U = exp[2iS/F n ] with 



i 



7T 



V 



V2 



T2 R ° " v¥? J 



(19) 



where F n is the pion decay constant. Under chiral rotations, U transforms like U — >■ LUR\ 
where L and R are SU(3) matrices belonging to SU(3)l and SU(3)r groups, respectively. 
To include the baryon octet, we introduce 



B 



J_vo , J_ A 
wo 



P 



n 



\ 



(20) 



7 



Again following the literature, we assume B transforms nonlinearly under chiral transfor- 
mation, 

B -> KBK ] (21) 
where K is a unitary matrix defined according to the transformation of £ = U 1 ^ 2 . 

£ -> L£^t, £ _> (22 ) 

It is clear that K is a nonlinear function of the Goldstone-boson fields. 

The quark-mass term breaks chiral symmetry and plays an important role in chiral ex- 
pansion. To exhibit its physical effect, the usual practice is to introduce the spurion field x, 
transforming as 

X -> L x R j • (23) 
However, to combine x with the baryon field B, we introduce x± 

X± = ± £x f £ , (24) 

which transform nonlinearly as x± ~^ Kx±K \ 

In the leading order, the chiral lagrangian for meson fields is 

C = -F%Tr[d ll rtd"U] + ]-F%BTr[M^U + U ] ' M] , (25) 
where M = diag{m u ,md,m s } is the mass matrix of light quarks. The leading-order chiral 



lagrangian for the baryon field is [24 



C = Tr ^Bi-fDpB - m BB + ^DB^J^, B + ^FB ltil5 [ Ufl , B] j , (26) 

where = i^d^ — is an axial vector current, D^B = d^B + [T M , B] and = 

{^d^ + 0n^}/2 is a vector current. 

A. Matching to CP-Odd Goldstone-Boson Operators 

Once there is a CP-odd term in the QCD lagrangian, it induces CP-odd terms in the 
effective Goldstone-boson lagrangian. These terms can annihilate odd-number (particularly, 
one) Goldstone bosons into the vacuum. Because of this CP-odd meson condensate, the 
original CP-even terms can now contribute to the CP-odd effects. Due to chiral symmetry, 
a meson condensate can generate physical effects only when the CP-even terms explicitly 
break the symmetry. 

As discussed in the last section, P-odd and CP-odd four-quark operators can be decom- 
posed into chiral (3,3), (6,6), (8,8) and their hermitian conjugate representations. They in 
turn can be matched to the corresponding chiral operators in the same representations. The 
leading ones without derivatives are unique and are shown in Table [TTJ 

We illustrate the matching process using O^f = uij^udd as an example. As discussed in 
the last section, this quark operator can be decomposed into irreducible representations of 
the chiral group using the spurion fields 

Q ud = ud,(3,3) + ud,(e,B) + ud,(8,8) + hc (27) 
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Rep. (3,3) (6,6) (8,8) 

Operator Of = iTr[Atf] Og 1 = iS%UfU) Of = iTi[H l UH 2 U^] 

TABLE II: Leading meson operators in individual irreducible chiral representations where A, S, 
Hi and H 2 are spurion fields in Table HI The appearance of i in front of each operator indicates 
that these operators generate P-odd and CP-odd vertices in the meson lagrangian; their Wilson 
coefficients in the lagrangian are defined to be real. 



Then, we can match each of the operators to the corresponding one in the meson sector 
through the non-perturbative Wilson coefficients C's 

ud,(S,3) ^ cC 3,3) 03 ) ud,m ^ cmo6 ^ ud,(8,8) ^ ^(8,8)^ (2g) 

The Wilson coefficients can be obtained by matching the simplest matrix elements: (0|O|7r°) 
and (0|O|?7), which can be calculated using non-perturbative methods such as lattice QCD. 

In this paper, we use factorization approximation to estimate these non-perturbative 
matrix elements. Lattice QCD calculations demonstrate that the matrix elements of four- 
quark operators can be factorized typically to within a factor of 2. Again take the operator 
as an example, which can annihilate 7r° and rj to the vacuum. [In principle, it also 
annihilates rj' , but this contribution is suppressed by the mass of rj'.] The annihilation 
amplitude can be estimated using vacuum saturation, 

(0\ui-f 5 udd\n°) « (0\dd\0) (0\ui-f 5 u\n°) . (29) 

Using chiral symmetry, one can get (0|uz7 5 m|7t ) = jr(0\uu\0) = —F^Bq, and (0|m75u|?y) = 

— FnBo/y/S. (This is consistent with the definition of the chiral rotation of U defined below 
Eq. (US).) Therefore, a term C A 0\{ in the QCD la grangian can be matched to the linear 
terms in 7r° and r\ in the chiral lagrangian 

C = C A BlFy + ^f A BlFl V + ... , (30) 

where ... represents higher-power meson fields. Then the leading terms in the potential of 
7r° and rj can be written as 



V = -B 
2 



[m u + m d )(7T ) 2 + ]-{m u + m d + Am s )rf 
o 



+ ^ mu _ md y v _ CiB 2 F 3^0 + ^_ r ^ j (31) 

which can be minimized to yield a condensate (7r°) and {rj). 

The above discussion can be easily generalized to an arbitrary four-quark operator, for 
which Eq. ( 130]) can be written as 

C = g n C,B 2 Fy + g v C A B 2 F^ + ... , (32) 

where g n and g v are numerical factors generated through the vacuum saturation approxima- 
tion. Then, the vevs of meson fields can be written as 

_ gpgCj [g^(m u + m d + 4m s ) - V3g ri (m u - m d )] 

A(m u m d + m d m s + m s m u ) 
, . _ BpF^Cj [-y/Sg n (m u - m d ) + 3g v (m u + m d )] 
N A(m u m d + m d m s + m s m u ) 
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which is inversely proportional to quark masses! The vev of U can be written as 



(U) = exp 







V 















2 

V3 



\ 



(v) J 



/F v 



(34) 



This defines the vacuum state of Goldstone-boson fields. 

Therefore, we can redefine the meson fields in the following way: 

U=(U)U' , 



(35) 



where U' collects the physical meson excitations. Through this redefinition, the meson 
lagrangian no longer contains terms annihilating the physical Goldstone bosons. Corre- 
spondingly, we redefine the baryon fields, 



(36) 



through a chiral transformation with L = (U) and R = 1. 

The above redefinition can change P-even and CP-even terms with explicit chiral symme- 
try breaking to P-odd and CP-odd terms. This is particularly true for the CP-even baryon 
lagrangian with linear dependence on quark masses, 



C c = dTr[BB]Tr[MlP] + c 2 Tr[M^BB^} + c 3 Tr[B^M^B] + h.c. 



(37) 



and 



C d = rf 1 Tr[5 75 5]Tr[M[/ t ] + d 2 Ti[M$} B lb B$}\ + d 3 Tr[B^M^B] + h.c. (38) 

Substituting (U) to the above equation, we get CP-odd pion-nucleon couplings through 

c{Vi[B' B'\Vi[{U) ] MU' ] ] + c 2 Ti[(U} j M^B'B'^} + c 3 Tr[B'^ (U) ] M^B'\ + h.c. (39) 

and the CP-odd masses of baryons 

diTr[B'j 5 B']Tr[(Uy ' MU' j ]+d 2 Tr[(U) ] B'j 5 B'^]+d 3 Tr[B'^ (U)^ B']+h.c. (40) 

which is part of the CP-odd mass generated by the four-quark operator. Note that since 
(U) is inversely proportional to the quark mass, the above contribution is not suppressed in 
the chiral limit. 

One can also get a CP-odd dipole moment by considering a photo-pion production term 
off the nucleon. When the pion is condensed through CP-odd effects, one generates a new 
contribution to the CP-odd moment, which is beyond the scope of this paper. 



B. Matching to CP-Odd Baryon Operators 

In this subsection, we construct the leading P-odd and CP-odd baryon operators induced 
by the CP-odd four-quark operators. These include all the operators with one baryon and one 
conjugate baryon fields, and without any quark masses or derivatives. All the independent 
operators are listed in Table IIHI 
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Rept. 



(3,3) 



(6,6) 



Operators 



of = iTr[BB]Tr[Atf],0. 



of = zTV[ J B 75 B]Tr[A?7 t ],0' 



W 

3 



(2) 
3 



o^isMmfiim, of 



= iTr[A^BB^},0 ( z } = 
iTi[A^B l5 B^],d^ 



iTr[Bf*A&B], 
= iTr[B^A^B] 



i^BBO^, 
Of = iSl](BOrm k m U l j7 Of = nv[BB]5gl7*^. 



(2) 



Q 
O 



of = *'&ES£tfr 1 ]' & [£ts£ff 2 ] j o. 



(1) = iTv[B^H 1 UH 2 ^B], Of 
(3) - iTr^BB^U^Hx], of 



)(4) 



=iTr[tB£H 2 gBgH 1 ], O, 



=iTi[B l5 ^H 1 UH 2 ^B], 6 



Of = iT^75.B£ff2l7tffi], Q 
Of =iTV[eS 75 e t ^i]Tr[^S^ 2 ], O 



(6) . 

(8) 
8 

(2) . 
8 

(4) . 



(6) 



Of =m[£BT 5 ZH 2 gBgH 1 ], 6 



n(8) 



iTt[B^H 2 WH^B], 

iTr[et J B 75 ^ 2 ]TY[ei?e t ^i] 



TABLE III: Hadronic operators that have the same quantum numbers as four-quark operators in 
different irreducible representations. 



There are two types of operators in Table IHIl those with and without tilde. For the 
first group without tilde, the expansion of the pion field generates the P-odd and CP-odd 
nucleon-pion vertices 



/-CP-odd 

NNtv 



(h c priTC + + h.c.) + h n nnn° + hpppn^ 



(41) 



For the second group, the leading order expansion is a bilinear-baryon term with a CP-odd 
mass structure, 

££as- s ° dd ~ -m*m l5 n. (42) 

This term contributes to the CP-odd baryon wave function. 

Traditionally, P-odd and CP-odd pion-nucleon couplings are defined in terms of isospin 



0, 1, and 2 of the operators, which can be written as [18 

Cnn = g^ N NT a Nn a + g^ N NN<f>° + g^ N (Nr a Nn a - 3Nr 3 Nn°) 



(43) 



where g^ NN is the coupling of the isospin-i term and r l are the Pauli matrices. Then, in 



terms of g^ NN , h c , h n , and h p can be written as 



h _ ^/2(o (0) +o (2) 
a c — v ^yg^NN-f-g^ 



\ \ ' • h n 



gnNNtg-n-NN'^gTrNNJ' lb P 



NN*g nNN 



as 



V&n) . («) 



where h p does not contribute to nEDM. 

To match the P-odd and CP-odd four-quark operators to the above baryon operators, one 
must find ways to calculate the corresponding non-perturbative Wilson coefficients. This 
can be done by considering the matrix elements of the quark operators in simple states. 
Take O^f = ui^udd as an example. As shown in the last section, it can be decomposed into 

group, 

+ u ^ m + 0^' {8,8) + h.c. 



irreducible representations of the chiral group, 



ud 



QUO 



o 



ud,(3,3) 



11 



(45) 
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The spurions related to this operator are given in Eqs. (JSJ, (|10|1 and (|T2|) . Q^ 3,3 ^ Q^i^fi) 
and (9"^ 8 ' 8 ^ must be matched to the hadronic operators in the same irreducible representa- 
tions and with the same spurions. Take the un-tilded hadronic operators as an example: 

3 

«J,(3,3) = £ C r(°OG ) +... ; 

i=l 
4 

^,(6,6) = ^ C W ( i ) + ... , 

i=l 
8 

^,(8,8) = ^ C W « + ... ; ( 46 ) 



i=l 



where represents higher order operators. 

Note that, an operator can be separated into hermitian part and anti-hermitian part. 
Since the QCD Lagrangian is hermitian, the hermitian part and the anti-hermitian part 
must have the same Wilson coefficient in the effective theory. Take the operator qiQR as 
an example, it is a (3, 3) operator, so it can be matched to CW in the chiral perturbation 
theory, while its hermitian conjugation qRQL is matched to CU with exactly the same Wilson 
coefficient since the QCD Lagrangian is invariant under the hermitian conjugate transfor- 
mation. Therefore, the hermitian part of qLQR can be matched to C(W + U)/2 whereas the 
anti-hermitian part can be matched to C(W — U)/2. As a result, one can use either the 
hermitian part or the anti-hermitian part of the operators to get their Wilson coefficients 
depending on which way is easier. For the operators without tilde listed in Table IHH the 
anti-hermitian parts contain terms having only one baryon field and one anti-baryon field 
which is easy to do the matching, while for the operators with tilde, the hermitian part is 
easier. Therefore, we choose to match the anti-hermitian part of the operators without a 
tilde whereas match the hermitian part of the operators with a tilde to get the Wilson co- 
efficients of them. One can show that this matching procedure works when current algebra 
is valid such as in non-relativistic quark model. 

The leading-order expansion of the hadronic operators are, for (3, 3) operators, 

{ 3 ] ~ ipp + inn + iAA + + + i£~£- + iE°E° + iE~E~ , 



Ojp ~ ipp + inn + — AA 



2i v 
3" 



Of ~ — AA + iE°E° + iE~E~ . (47) 
3 

Therefore, we can determine the Wilson coefficients with four physical matrix elements, 

qw + cf = (-z)(p|orf (3 ' 3 V> , 

= (-z)(S°|0 1 f (3 ' 3) |S°) , 
C^ + fcf = (-.)(A|0^ 3) |A) , 
tfW + Cf = (-t)(E°\0^' m \E°) , (48) 
where we have chosen the normalization condition 

(P\P') = (2tt) 3 5 3 (P - P') , (49) 
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where P and P' are the momenta of the states. 

Since the number of equations is larger than the number of variables, to get a solution 
the following condition must be satisfied, 

/no (p\o^' 3) \ p ) \ 



det 



1 (S°|On' {3 ' 3) |S°) 
1 | (A\0^' m \A) 
V 1 1 (S°|On' (S ' 3) |H°) / 







(50) 



which gives a nontrivial relation among these matrix elements; 

_ -( p \O ud ^\p) - -(E |O"f' {S,3) |E > + (A|Oj? (S,3) |A) 
3 3 



. 



(51) 



This relation must be satisfied in the chiral limit, so it is a test for direct calculations of the 
matrix elements. Similarly, a simple inspection of Eq. (I47p can give us some more relations 
among matrix elements 



(p\0^' 3) \p) = (n\0^' m \n) 

0|r>«d.(3,3)|£,o\ 



(E u |0 
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(E°\0^ m \E°) = (E-\O l u 

Generalizing the above discussion to (6, 6) and (8, 8) operators, we write down the leading 
expansion of the hadronic operators, 



(S + |On' (3 ' 3) |S+) 



(£-|On' (3 ' 3) |£-> , 



(52) 



~ -AA-*E°E 
3 



(2) 

o£ 3) - -AA + /E"E 



iE + E+ + iE"E" 



zpp + inn + -AA + zE°E° + iE + E+ + zE^E" 



z - 



zE°E° + iE + E + + iE"E" 



2i (kp + nn + AA + iS°S° + S~S~ + E°E° + E+E+ + E~E~) 



(53) 



^8 

^8 

(3 
•^8 

(4 



o, 
o 



(r 

w 8 



•^8 



/ I -AA 

6 



' AE + ^E A+is E + S + E- 
2VS 2 



1 



1 



O'-'O 



/ I -AA + EE 



2^3 
1 



AE C 



AE U 



i(nn + -AA 3_AE° 

iE+E+ , 



^E A + -E E° + E-E- 
2a/3 2 

_^E°A + iE E° + E + E + 
2^ 2 

-^=E°A + -E E° + S-E- 
2v/3 2 



~ i I -AA + -^=AE° - E°A - -E U E 



1. 



2^ 



/ ( iAA - ^AE° + E°A - -E°E° 
6 2%/3 2 



(54) 
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from which we can get similar relations among matrix elements just like in the (3, 3) case 
shown in Table IIVI The other four-quark operators with the same flavor structures have the 
same relations among hadronic matrix elements as in this case. 



Rep. 


Relations 


(3,3) 


(p\O u ^ (3 ' 3) \p) = (n|O^' (3 ' 3 V),(S |On ,(3 ' 3) |H°) = (H^f (3 ' 3) |~->, 
(S°|0^' (3 ' 3) |S°) = (S+|0^' (3 ' 3) |S+) = (S-|0^' (3 ' 3) |S-), 

2 /„J/n urf >( 3 > 3 ) U\ 1 /v0l/o ud >( 3 > 3 ) I v0\ i / A l/O urf '( 3 ' 3 ) 1 A \ n 

-a\P\°n \P)-3^\ ii I s ) + \ A \°ii l A > =0 


(6,6) 


{p\C$<fi>%) = (n|0^' (6 ' 6) |n),(S + |On' (6 ' 6) |S+) = (E~\0^' (6 ' 6) |E"), 

(h + |o w ^ 6 ' 6 ^|h + ) = (s~|o"^ 6 ' 6 ^|e; - ) 

2{p\O u 1 f'%) - (S°|0^' {6 ' B)) |S°) - 3(A|0^' (6 ' e) |A) +2(H |O^' (6 ' B) |S°) = 


(8,8) 


{p\0^' {8 ' 8) \p) + {n\0^ 8) \n) + (S°|0^' (8 ' 8) |S°) - 3(A|0^' (8 ' 8) |A) 

+ <~ |O 1 f' (8 ' 8) |H°> + (E-\0^' {8 ' 8) \E-) = 
(p\0^ {8 ' 8) \p) - (n\0^' {8 ' 8) \n) - (H°|0^' (8 ' 8) |H°) + (Z~ \0^ {8 > 8) \~- ) 
-V3(A\O^ 8) \^°)-V3(J: \O^ {8 ' 8) \A) = 



TABLE IV: Relations among hadronic matrix elements of the four-quark operators in different 
chiral representations. 

One can either build models or do lattice QCD calculations to get these simplest four- 
quark matrix elements. Once known, one can get the Wilson coefficients by solving Eq. ( 148]) 
and similar equations for (6, 6) and (8, 8) operators. Then one can expand these hadronic 
operators to the first order with one meson field in each term to get the P-odd and CP-odd 
pion-nucleon vertices. A similar method works for baryon operators with tilde. We will 
consider these matrix elements in the next section. 



C. Matching to EDM-Type Operators 

In ;\PT, any CP-odd quark-gluon operator will generate directly an EDM contribution 
to the neutron, analytical in the chiral limit. To write down such a contribution, introduce 
vector and axial vector octet potential and a M , which transform under local chiral rotations 
(with space-time dependent chiral transformation) as 

J> = + a M — > v'^ + a'^ = R(v^ + a fJ )R j + iRd^R ] , 
In = — > v'^-a'p = L(v^ - a M )L f + iLd^l) . (55) 
The corresponding gauge fields are defined as 

/£, = d^-dJ^-i^Q . (56) 
The gauge fields with definite parity are defined as 

f% = ?fZt ± Uiut\ (57) 
which transform under chiral transformation as 

ft Kf±K* . (58) 
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When reducing to the electromagnetic field, a M = 0, f^ v = (&Q£ ± ^,Q^)F^ U , where Q = 
diag(2/3, —1/3, —1/3) and F^ v is the electromagnetic field [24 . 

One can write down a number of EDM type of operators which contain B and B, f^ v , and 
the spurion fields A, H, and S. These contributions are direct matching contributions to 
the neutron EDM, and cannot be calculated in %PT. These chiral constants can in principle 
be calculated in lattice QCD. However, we will present quark-model estimates in Sec. V. 



D. Peccei-Quinn Symmetry and Induced 0-Term 



The experimental upper bound on the neutron EDM gives a strong constraint on the 

Using the current 

(59) 



P-odd and CP-odd 0-tenn, 6GG, in the QCD lagrangian [TJ, [l4j, [15 
experimental limit [H^ 



d n < 2.9 x 10 



-26 



e cm 



one can get the upper bound, 



6 < 10" 10 



(60) 



On the other hand, it is unnatural for a parameter of the fundamental theory to be so small 
without fine tuning. There are generally two ways to solve this strong CP problem in the 
literature. The first is by introducing the spontaneous breaking of parity. Since the #-term 
also breaks parity, if at some high energy scale parity is conserved, then the #-term at low 
energy scale can only be generated by loop effects and will be suppressed naturally (3lj . 

The other way is to introduce the Peccei-Quinn symmetry, U(1)a (32[. After the sponta- 
neous breaking of the symmetry, there emerges a pseudo-goldstone boson, a, which is called 
the axion 33M35| . The effective Lagrangian for the axion field can be written as 



2 / a 87T ^ 



(61) 



which includes an effective interaction with GG. The axion field gets a small mass through 
the non- vanishing correlation function 



K = i 



d A xe ik - x 











k=0 



(62) 



after taking into account the non-perturbative QCD effect (13, 35| . 

When there is an additional neutral P-odd, CP-odd quark operator, Ocp-odd, in the 
lagrangian, the correlation function 



d 4 xe ik - x ( 



T^GG(x), Ocp-odd(O) 







(63) 



fc=0 



will be generally nonzero. Therefore, the vev of a, which cancels precisely the #-term in 
the original lagrangian, will now be shifted by a small amount proportional to K\. A non- 
vanishing effective #-term is induced as 36 



hid 



Ki 
K 



(64) 



which can contribute to the neutron EDM. 
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Following Ref. 36], we take the operator ui^udd as an example to calculate the con- 
tribution to neutron EDM through the induced 0-term. Then, K\ can then be written 

as 



K x = % 



d 4 xe ik - x 







T ( — GG(x) } C^ui^uddUS) 
8n 







k=0 



Using the chiral anomaly [37J, one can get 



47T 



GG = d^J§ — 2m*(m75M + dij$d + si^f^s) 



where 



Then one can get 



J' 1 



m 



-M7 M 7 5 w H -c?7 M 7 5 d H -37^753 

m u m d m s 



(65) 

(66) 
(67) 



aW*-* <0 |T (d^{x), C 4 ui-f 5 udd(0)) I 0) 



fc=0 



d A xe ik ' x (0 |T (2m t (w7 5 u + c% 5 <i + 34753) (ar), C 4 ui>y 5 udd(0)) | 0) fc=0 . (68) 



The second term on the right-hand side of the above equation is negligible compared to the 
first term because it is explicitly proportional to the reduced quark mass m* and the operator 
uij 5 u + di^d + sz7 5 s cannot annihilate light mesons. Therefore K\ can be calculated as 



K x « l - J d*xe tk - x (0\T(d^(x),C 4 m l5 udd(0))\0) k=0 



--C 4 (0\[Q 5 (0),ui l5 udd(0)}\0), 



(69) 



where Q5 i s the charge related to the current Ji* defined in Eq. (1671) . In the spirit of large 
Nc Jlsl 39] expansion one can assume that 



(0\uij 5 udij 5 d\0) < (0\uudd\0) « <0|uu|0)<0|dd|0>. 



(70) 



Therefore, we can get 



7 — 171 — 

K x ~ — C 4 (0|[Q 5 (0),m7 5 M(0)]|0)(0|dd|0) = -C 4 — (0|uu|0)(0|dd|0) 
2 m u 



m,. 



Using the previously known result [35 



A' 



one can get the 9 angle induced by the operator uij 5 udd, 



find — TT 



C^BqF 



K m u 

A similar result can be obtained for any other CP-odd four-quark operator. 



(71) 



(72) 



(73) 
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Using the standard chiral result in the literature (13|, we write down the effective chiral 
lagrangian corresponding to this induced 9 term; 



A0m - - - 

C e = -=-^ (c 2 Tr[EBB\ + c 3 Tr[BEB\) + 2m,e{3d 1 + d 2 + d 3 )Tr[Bi l5 B] . (74) 

From the above, we read off the CP-odd pion-nucleon coupling and the CP-odd mass of the 
neutron; 

2y/2CjB F v m* _ 2C±B F 7T m* 

fin fin. i 



M t = 2C * BoF * m * ( 3dl + d2 + d3 ) . (75) 



m,. 



Comparing this with the meson condensates contribution in Eq. ( 177|) . one finds that they are 
in the same order. If the Peccei-Quinn symmetry exists, one should add this contribution 
to the neutron EDM. However, since it is not known if the axion mechanism is in operation, 
we will not include this contribution to the nEDM in the remainder of the paper. 



IV. P-ODD AND CP-ODD NUCLEON-PION VERTICES AND CP-ODD NU- 
CLEON MASS 

In this section, we study the induced physical P-odd and CP-odd nucleon-pion vertices 
as well as the CP-odd nucleon mass from four-quark operators. There are a number of con- 
tributions to consider: First, the CP-odd meson lagrangian will generate meson condensates 
which can convert a CP-even vertex into a CP-odd one. Second, the baryon wave function 
contains the CP-odd part due to the CP-odd nucleon mass, which can also rotate a CP-even 
coupling into a CP-odd one. Finally, there is the contribution from the direct matching 
operators (without a tilde) in TABLE II. We will consider all of these in this section. 



Meson Condensates Contribution 



We use the vacuum saturation approximation to calculate the meson effective lagrangian; 
the vevs of 7r° and 77 can be obtained from Eq. (1331) . where g n and g n for all the four-quark 
operators built with color-singlet and octet scalar currents are listed in Table |V] Those 
induced by tensor operators vanish in this approximation. 



In the large N c QCD [38[ (also see Ref. [39| for a good review), the leading contributions 



for operators constructed from two color-octet currents and two tensor currents are shown 
as diagrams (a) and (b) in Fig. [TJ respectively. Detailed analysis shows that the diagrams 
(a) and (b) suffer from 1/N% suppressions compared with (c), which stands for the operator 
constructed from two scalar color-singlet currents. 

Terms contributing to the P-odd, CP-odd nucleon-pion vertices through the condensates 
of neutral mesons are shown in Eq. fl3"§l) . At tree level, one can relate the coefficients Cj, c 2 , 
and C3 to the mass differences of the baryons and the ttN a-term, and their values can be 
found in the literature 124 



Cl = 2B b , c 2 = 2B (b d -b f ) , c 3 = 2B (b d + b f ) , (76) 
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Operator 


9-k 


9r) 


Operator 


9-K 




ui^/^udd 


i 


1 1. /q 


ui^j^t a udt a d 


U 


U 


di'j^duu 


— i 


1 / /q 


di^§t a dut a u 


U 


n 
U 


ui^uss 


i 
i 


1 /. /q 
1/Vo 


ui^fc ) t a ust a s 


U 


U 




n 
U 


/. /q 


si r y^t a sut a u 


U 


U 


di^ftydss 


— 1 


1 /-/q 


GL%^§1 QjSL S 


n 
u 


n 
u 


si'jssdd 





-2/x/3 


si^t a sdt a d 








ui'Ysuuu 


5/6 


5/(6 V3) 


ui^^t a uut a u 


-2/9 


-2/(9^3) 


di^/^ddd 


-5/6 


5/(6^3) 


di^t a ddt a d 


2/9 


-2/(9>/3) 


sij^sss 





-5/(3 V3) 


sij5t a sst a s 





4/(9^3) 



TABLE V: ^ and induced by four-quark operators constructed by scalar currents. Those 
induced by products of tensor currents are zero. 





(a) 



(b) 




(c) 

FIG. 1: Annihilation of pion by four-quark operators: (a) operator constructed from two color- 
octet current, like ui^t a udt a d; (b) operator from two tensor currents, like ua^i^f^uda^d; (c) 
operator from two scalar currents, like uij^udd. 



where b = -0.517 GeV -1 , b d = 0.066 GeV" 1 and b f = -0.213 GeV" 1 . 

The vertices we are interested in have two nucleons and one pion because of the infrared 
enhancement in the pion loop [15]. From Eq. ( 139]) we can read off the relevant terms, 

— {c 3 [3v / 2(vr°)(m u - m d ) + V6(r])(m u + m d )\ J (npTr~ + pnix + ) 

, ;/ , 2 ^c 3 m d (3(TT°) - y/3(ij)) + Ci[3(m u + m d )(rr ) + V3(r))(m u - m d )]| nmr° , (77) 



3F3 
2 



in which (71*°) and (77) are given in Eq. ( 1331b It is customary to define the P-odd, CP-odd 
nucleon-pion couplings 



£-cpv — h c (pmr + + npir ) + h n nmi , 



(78) 



where h c and h n induced by meson condensates are listed in Table IVTl Typical values of h c are 
one order of magnitude smaller than the value of h n because \/2c 3 (m d —m u ) <C Aci(m u +m d ). 
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TABLE VI: CP-odd pion-nucleon couplings induced by meson condensates. C4 is the Wilson 
coefficient of the corresponding four-quark operator. The two columns on the right side shows the 
P-odd and CP-odd pion-nucleon vertices calculated using factorization which will be discussed in 
Sec. V. 

For h c or h n generated by a certain four-quark operator, if the contribution from (tv°) is non- 
vanishing, the contribution from (rj) can be neglected since (tt°) / (77) ~ m s /m ~ 30. This also 
explains the contributions from operators with the sij^s factor are much smaller than those 
without. Finally, the contributions from operators made of color-octet currents are smaller 
than those from operators made of color-singlet currents because a Fierz transformation is 
needed for color-octet operators to annihilate the mesons, introducing a suppressing factor 
of 1/4. 

In Table I VI I one can see that the P-odd and CP-odd pion-nucleon couplings are propor- 
tional to Bq, which is related to the quark condensates. The value of Bq can be extracted 
from the pion mass 

ml = B (m u + m d ) . (79) 
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The natural scale for x?T is |40j , and for simplicity we use the same scale to define the 
quark masses to get Bq. The quark masses we use are m u = 2.4 MeV and = 4.75 MeV 
in MS at 2 GeV. Using the one- loop renormalization group to run them down to fi = AttF^, 
we have 

Bo = 2.2 GeV . (80) 
Here we have used one-loop Aqcd = 250 MeV. 



B. Direct Contribution from Matching 



To get the P-odd and CP-odd meson-nucleon coupling through direct matching, one 
needs to calculate the matrix elements listed in Table IIVI Lattice QCD is perhaps the 
ultimate choice for calculating hadronic matrix elements. However, it is still quite difficult to 
directly calculate the matrix elements of four-quark operators between baryons. Therefore, 
we resort to quark models to get an estimate. In the remainder of this subsection we will 
use two different quark models to calculate these hadronic matrix elements: the simple 
non-relativistic quark model 41 - 45 1 and the MIT bag model 45-4^|. We also discuss the 
significance of the model calculations from the viewpoint of naive factorization. 



1. Non-relativistic Quark Model 

Here we consider the simplest version of the non-relativistic quark model with harmonic 
oscillator interacting potentials, 

h = -J2^ + 1 ! y uj2 " f2)2 + (rl " fs)2 + (fs " fl)2] ' (81) 

where fx, f 2 , and f 3 are positions of the three quarks inside the baryon, m c is the mass of 
the constituent quarks and u is the angular frequency. One can isolate the center of mass 
by introducing the Jacobi coordinates, 

1 

R = -^(fi + f 2 + f 3 ) , 

P = -j=(fi-r 2 ) , 

A = _(ri + 7*3-2^3) • (82) 

Then the spatial wave function of the nucleon can be written as 

/(f 1; r 2 , r 3 ; k) = (3V3)- 1 / 2 exp(zP ■ R/V3)ip(p, A) , (83) 



where if?(p, A) = (a 3 /ir 3 / 2 ) exp [-a 2 (p 2 + A 2 )/2] in which a = (mu) 1 / 2 w 0.41 GeV |4| is 
the oscillator parameter, and P is the nucleon momentum. It is easy to check that the 
wave function is normalized to (27r) 3 <5 3 (P — P'). The internal part of the wave function is 
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assumed to have SU(6) spin-flavor symmetry. For example, the spin-up proton state has 
the following wave function, 



|Pt> 



d 3 nd 3 r 2 d 3 r 3 f(f h r 2 ,r 3 )e abc \uf (n)d\\f 2 ) - (n)d b ^f 2 )] v*(f 3 )\0)(M) 



where a, b, and c are color indices and the anti-commutation relation of the non-relativistic 
quark creation and annihilation operators is defined as {u^ (x) , u b Jy)} = S a b5 a p8 3 (x — y) 
with a and as spin indices. The spatial part of the wave functions is common for all 
members of the baryon octet. The SU(6) internal wave functions are listed in Table rvTTI for 
easy reference. 
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TABLE VII: SU(6) wave functions of baryon spin-1/2 octet. 

Using Eqs. (I7|)- (1T21) . one can project operator ui^udd into different irreducible represen- 
tations of the chiral group, 0^' i3 > 3 \ 0^' (6 ' e) , and 0^ (8 ' 8) as in Eq. (Hoi) . Restricting to th< 
non-relativistic case, these operators become 

^m {x) „ _i_ . ^ x)<{x) ^ {x) ^ x) _ d a J(x)u a a (x)u^(x)d b (x) 
'u a J(x)u a a (x)d b ^(x)d b 3 (x) + df (xK(x)«J(x)4 (*) 

l /3 \^! a j3 tn ' * 



0^ (6 ' 8) (x) 



^11 



Md,(8,8) 



~ - : <t (x)<(x) ^t (x) ^ (x) . 



(85) 



where w and c? are non-relativistic two-component quark annihilation operators, a and h label 
the color, a and /3 label the spin, and the ": :" means that the products of the constituent 
quark fields are normal-ordered. 

Considering the (6, 6) component as an example, the simple quark model gives the fol- 
lowing matrix elements: 



(j>,{p)\ot' m \ P ,{P)) 



(n,{P)\Ot' m \n,{P)) 



(s+(P)|o^ i(6 ' §) |s+(P)) = (£-(p)|o^ (6 ' S) |£-(p)) = o 

(E (p) | O ^(6,6) |E (p)) = 

(At(p)|0 ^(6,6) |At(p)) = Q) 

(s (p)|O ^(6,6) |H (p)) = (St - (P) |0^(^)| 5f(p)) = 



(86) 
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where a = J d 3 r f* (P; x, x, r) f (P; x, x, r) is independent of x. It is easy to check that these 
matrix elements satisfy the symmetry conditions listed in Table HVl Using Eq. (I53p . one 
can get the Wilson coefficients for (6,6) hadronic operators defined in Eq. (j4"6"l) ; 



r (i) _ r (2) _ 1 a 3 



8(2tt) 3 / 2 ' 

cf ) = ef = . (87) 

Expanding the hadronic operators to the first order, one can get the P-odd, CP-odd three- 
point nucleon-pion couplings, h c and h n . The result induced by O^ 6,6 ) i s 

h c = C^c?/(%^ /2 Fk) ~ 0.022C , r i V/P 7r , K = ■ (88) 

In the same way one can calculate h c and h n induced by the (3, 3) and (8, 8) components 
of ui'jsudd. Taking into account the hermitian conjugate part of each component, the 
contribution for h c and h n is doubled. 



2. MIT Bag Model 

The basic idea of the bag model is that valence quarks are confined in a bag where the 
vacuum is in a phase different from the true QCD vacuum. The inside has a constant 
energy-momentum density generating a negative pressure, B, which is balanced by the 
positive pressure of the quarks. The bag is usually taken as a sphere of radius Rq. The 
quarks inside the bag move freely with the following wave functions, 



fc,_ 1 , 1/2 ,„(f,i) = ^( «<^W«°)x ? , 

V4vr V -Ji(u n ,-ir/Ro)(T • rx r , 



The normalization factor of the above is 

,3 \ V2 



N M = ^37 T\—2 ) ' (9°) 

\2R*(u nK + k) sin ui nK J 
The boundary condition gives the energy eigenvalue equation, 

tanw nK = UnK , (91) 

UJ nK + K 

and numerical calculation gives Uq = 2.043. The ground state of quarks is K — — 1, n — 
state. For the baryon octet, all the quarks are in this state. Keeping only this, the quark 
operator can be written as 

q(x) = ipo,-i,i/2,rn(x)e~ iuJO '- lt/Ro b ^ 1A/2 , m + (anti-quark creation) . (92) 

The physical meaning of the operator 6 m (0) is that it annihilates a quark with quantum 
number described by the wave function i[)a-i,i/2,m- Due to the assumption that inside the 
bag the interaction between quarks and gluons is negligible, flavor and spin automatically 
become good quantum numbers. 
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TABLE VIII: P-odd, CP-odd three-point pion-nucleon vertices generated by P-odd, CP-odd four- 
quark operators. The couplings induced by operators constructed by two color-octet currents are 
equal to the couplings induced by corresponding color-singlet operators multiplying by —2/3. 



We again take the (6, 6) component of uij^udd as an example, which can be written as 



+ d a Md\ 



(93) 



where a and h are indices of color, A, <r, p, r labeling the spin. The creation and annihilation 
operators here are just like &o,-i,i/2,m m Eq. f )92|) . Using 



N 2 



57T 



K(W7-Ro) - 3i(uor/R )] S X a - 2ij (ujr/Ro)jo(ur/R)xxV • ^Xo 



(94) 



and only keeping the terms which give non-vanishing contributions after integrating over a 
spherical region, we have 



N 4 

64tt 2 v 

-Aj^(u r/R )j 2 (uj r/R )(x{a ■ rx*)(xl* • rXr] 



[j%(u o r/R ) - j 2 {u r/R )] 2 5 Xa 5 pr 



where we neglect the term proportional to a ■ f. In a proton state normalized to our con- 
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vention before, the expectation value of the operator can be written as 

(Ptl< 6,e) bt) = -^K) 4 ^/^ 



x |[j 2 ( Wo r/i?o) - JxWARo)] (Ptl^bt) -±ti{^r/Rv)ji{uvr/Rv)(^\B\^)^ , (95) 

where 

A = : u a Md b Jd b p : + : d?«dj : 

B = : <t (x t«7 • r Xc )uld b }{ X \& ■ f X r)d b T : 

+ : d a \ X \ B ■ rx„)«(xj<? ■ f Xr )< : • (96) 

A straightforward calculation gives 

(p t \A\p t ) = (n t |A|n t ) = 1 
(E+|A|E+) = (Ef |A|Ef) = 

(E?|i4|Ef> = 2 
(A t |A|A t ) = 

(H°|A|S°> = (Sf |A|8f> = , (97) 



and 



(pt|B|pt> = (n t \B\n t ) = 1/3 
= (Ef|S|Ef) = 
(E^IE^) = 2/3 
(A t |5|A t ) = 

<S°|B|S°> = (HjT|S|Hf> = . (98) 



Therefore we can get in the MIT bag model 



( Pt |Orf 6 ' §) (x)b t ) = *A + ^B, (99) 



where 



" = \ NM 'h^l (lk) d (i^) ^("or/BoJjf^or/flo) , (100) 

and similarly for other matrix elements. 

Then, using the method we used in the non-relativistic quark model, we can get h c and 
h n induced by O 1 ^ , 

^ = ^^' kn = °- (101) 

One can compare this with the result from the non-relativistic quark model in Eq. (|88|) . 
where h c is proportional to a 3 . From the definition of ip below Eq. (153"]). 1/a can also be 
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TABLE IX: Same as Table lyTPIj except the matrix elements are quoted here at the scale ^ = ^ttF^ 
assuming the quark model scale of 400 MeV. 

seen as the radius of the baryon. It is well known that 1/a = 0.5 fm gives a too small value 
for the proton's charge radius and the pion cloud is usually invoked to gap it. On the other 
hand, the bag radius is usually taken to be 1.0 fm, which will give a considerably smaller 
h c . In any case, it is reasonable to consider Rq ~ 1/a and take the non-relativistic quark 
model result as the representative. 

The couplings h c and h n induced by color-singlet four-quark operators are listed in Table 
IVIIII and those by color-octet operators are equal to the above multiplying by —2/3. In 
Table IVHH many four-quark operators yield zero h c and h n because we neglect the "sea 
quark" contribution. By making the four-quark operators normal ordered in Eq. (185]) and 
(|96|) . one cannot get any contribution to h c and h n from four-quark operators containing 
strange quarks. 

Model calculations do not have explicit QCD scale dependence. To match the results 
with QCD matrix elements, we have to assume a model scale and using perturbative QCD 
(pQCD) evolution to run them to appropriate perturbative scale, for which we choose to be 
fi = An Fx- In this work, we assume the model scale to be at 400 MeV and Aqcd — 250 
MeV and take into account the pQCD effect using one-loop renormalization group equation 
to run the operators down to the energy scale of the model. At this low energy regime the 
strong coupling is large and the one-loop pQCD evolution is by no means accurate, but it 
may still serve as an estimate of the pQCD effect. The matrix elements at scale \x are shown 
in Table ED 

3. Contribution from odd-parity resonances 

The P-odd and CP-odd quark operators can also generate a CP-odd pion-nucleon in- 
teraction through the parity-odd excited resonances which is shown in Fig. [2j The P- 
odd and CP-odd quark operators can generate mixings between nucleons and parity-odd 
excited resonances which can be calculated using quark models [H(|. Take the operator 
Oyf* = ui^^udd and the intermediate state N(1535) as an example, using the harmonic 
oscillator non-relativistic quark model the mass mixing between neutron and A^(1535) res- 
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7T 



N N*,A* N 

FIG. 2: P-odd and CP-odd pion-nucleon coupling generated by the four-quark operators through 
parity-odd resonances, where the black dot is the CP-odd, four-quark operator, N* and A* are the 
CP-odd excited states. 



onance can be estimated as 8 = m c oo 2 /(8 v / 3vr 3 / 2 ), where m c to w 300 MeV are the 
constituent quark mass and the frequency of the harmonic oscillator, respectively. The res- 



onance can decay into a nucleon plus a pion, the partial decay width is about 50 MeV |51 
The effective Lagrangian for this process can be written as 

C N , = g N *NN*jc + h.c. , (102) 

where as an order-of-magnitude estimate we discard the isospin quantum number. Then, 
from the partial decay width one can get g^* ~ 0{1). The P-odd and CP-odd pion-nucleon 
coupling induced by this mixing can be written as 

h mix = ^ 49N * 6 « 6 x 10- 4 C 4 GeV 2 , (103) 
M N » - m n 

where C4 is the Wilson coefficient of the four-quark operator. Compared with the direct 
matching contribution listed in Table IVIIH one can see that h m i X is about two orders of 
magnitude smaller and therefore its contribution to nEDM is negligible. 

The contribution from Fig. [2]can be seen as a one-loop contribution since the intermediate 
resonances may also be described as scattering states of pion and nucleon. Therefore, this 
contribution is suppressed by a loop factor. 



C. Tree-Level CP-Odd Mass of Neutron 

The nucleon CP-odd observables receive contributions from its CP-odd mass term 
ml 'fajsif) . In %PT, there are also two sources of CP-odd mass: that induced by the conden- 
sates of meson fields, namely (ir°) and (77), and that from the direct matching contribution 
of the four-quark operators. 



1. Meson Condensates 

The relevant terms contributing to the CP-odd mass of neutron can be read from ex- 
panding Eq. (jlDl) . which gives 

If 1 2 

ni~f 5 n— l-di[(m u - m d )(iT ) + y|( m « + m d)(v) ~ -^ m s(v)} 

+d 2 ^=m s ( V ) + d 3 m d ((n°) - , (104) 
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where di, d 2 and d 3 can be related to the discrepancy of the Goldberger-Treiman relation, 
and the values d<i and d^ have been determined in the literature [521 . 



d 2 = -2B m (D l9 - F w ) , d 3 = -2B Q m (D w + F l9 ) 
where mo is the common octet mass in the chiral limit, and 

m F 19 « -0.2 , m D 19 « -0.4 . 



(105) 



(106) 



Note that the signs of the _F 19 and D± 9 here are different from those in Ref . [52| . Since d\ has 
not been determined from isospin- violation effect, we will set it to be zero in the following 
calculation. One should note that disregarding d\ leads to some errors because m s (rj) might 
be the same order as rridi^ )- 



Operators 


m' n /(W- 3 C 4 B$ GeV) 


Operators 


m'J (KT 3 C A Bl GeV) 


ui^udd 


-8.8 


ut a i^udt a d 





di^^duu 


5.7 


dt a i'y 5 dut a u 





ui^dss 


-8.8 


ut a i'y^ust a s 







3.2 


st a i^§sut a u 





di'jsdss 


5.7 


dt a i^dst a s 





sij^sdd 


3.2 


st a i^5sdt a d 





uij^uuu 


-7.4 


ut a i r )z > uut a u 


2.0 


di^^ddd 


4.7 


dt a ij 5 ddt a d 


-1.3 


si^sss 


2.6 


st a ij5sst a s 


0.7 



TABLE X: CP-odd mass of the neutron induced by meson condensates. Contributions from oper- 
ators made of tensor currents are neglected due to the large- Nc suppression. 



2. Direct Contribution 

The leading-order expansion of the tilded hadronic operators listed in Table II III are her- 

~ (2) 

mitian. Take & as an example. It can be written as 

Of ~ pi l5 p + ni lb n + + S i 75 S + S+z 75 S + + fri^TT , (107) 

which gives a CP-odd mass of neutron. To calculate the matching coefficients, we can see 
from above that the leading-order expansion is parity-odd, and we need to calculate a parity- 
odd quantity. The simplest is As ■ Ap, where As* is the spin difference between the initial 
and final states and Ap is the momentum difference between the initial and final states. 

In the non-relativistic quark model, take the (6, 6) components of ui'j^udd, as an example, 
to calculate the matrix elements proportional to As- Ap; the relevant part of the four-quark 
operator can be written as 

Ot{ m - -\ : [V • {u^u)\{d}d) : +£- : (^)[ V • (Sad)} : 



' : [V ■ (d^au)](uU) : : (rf f n)[V ■ : \ , (108) 



2m c ' 2m c 
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TABLE XI: CP-odd mass of neutron induced directly by color-singlet four-quark operators. The 
CP-odd mass induced by color-octet four-quark operators are equal to the one induced by corre- 
sponding color-singlet operators multiplied by —2/3. 



Operators 


CP-odd mass/(a 3 C4) 


Operators 


CP-odd mass /(aPC^) 


ui'jsudd 


0.212 


ui^a^uda^yd 


0.0280 


di^duu 


-0.336 


ui-f 5 a^t a uda fiv t a d 


0.189 


uij5t a udt a d 


-0.0314 






di^t a dut a u 


0.0799 






ui^uuu 





ui^5t a uut a u 





di^^ddd 


-0.249 


di-/ 5 t a ddt a d 


0.0968 



TABLE XII: Same as Fig. IXII The matrix elements are now evolved to the scale where \x = 4nF n . 

where u and d are two-component quark operators, mp is the mass of the constituent quark 
which is set to be one-third of the nucleon mass. The wave functions of baryons in the non- 
relativistic quark model are listed in Eq. f l84"|) and Table IVHI Then using the same method 
as described in the last section one can get the CP-odd mass of the neutron directly induced 
by the tilded operators, and the results are listed in Table IXII After the leading-order QCD 
evolution to the scale where \i = 4:irF n , the result is shown in Table IXlTI 

3. Contribution to CP- Odd Meson-Nucleon Coupling 

If rotating away the CP-odd nucleon mass through ?7a(1) transformation, one can gener- 
ate new contributions to the CP-odd meson-nucleon coupling from CP-even chiral operators. 
However, this contribution is of higher order in chiral power counting because all the CP-even 
meson-nucleon interactions are suppressed in the chiral limit, whereas the CP-odd coupling 
we considered in the previous subsections are not. 

V. FOUR-QUARK CONTRIBUTION TO NEUTRON EDM IN X PT 

In this section, we study the CP-odd four-quark contributions to the neutron EDM in 
%PT. The approach here is completely general and is applicable to any CP-odd quark-gluon 
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operators. Some results presented can be found in the literature; however, to our knowledge, 
this is the most systematic and thorough discussion in the context of the CP-odd four-quark 
operators. In the last subsection, we make a comparison of the four-quark contributions in 
different approximations of non-perturbative QCD physics. 

In xPT, the leading contributions come from many different sources. Since the CP- 
violating pion-nucleon couplings are 0(1), the pion loop contribution to the neutron EDM 
is 0(1), apart from possible enhancement by chiral logarithms. On the other hand, the 
direct matching contribution is also 0(1), along with the pion condensate contribution 
through photo-production amplitudes. Finally, the CP-odd mass terms contribute through 
the nucleon magnetic moment after chiral rotation. This contribution is again 0(1) in chiral 
power counting. We will consider all these leading contributions in the following subsections. 
We ignore the subleading contribution in this work. 



A. Direct Matching from Quark Model 



We have first considered the direct matching contribution from the four-quark operators 
to the neutron EDM in Sec. IV. When any CP-odd quark-gluon operator is matched in xPT, 
there appear many tree-level neutron EDM-like operators in the chiral Lagrangian [ 1 3| . We 
do not have much to say about the size of the Wilson coefficients other than they are 0(1) 
in chiral power counting. Since they also serve as the counter terms for ultraviolet-divergent 
chiral-loop calculations, they depend on the regularization scheme and subtraction scale. 
In this work, we choose to estimate this contribution using nucleon models with dipole 



excitations into odd-parity resonances, such as Su, following the work in [50 



N ,A 



FIG. 3: Direct calculation of the neutron EDM in quark models. The neutron makes a transition 
to a CP-odd excited state and goes back via electromagnetic interaction, where the black dot is 
the CP-odd, four-quark operator, N* and A* are the CP-odd excited states. 



We use the non-relativistic quark model with harmonic oscillator potentials to estimate 
the contribution from the first CP-odd excited states, which is shown in Fig. [3j The wave 
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functions of the lowest CP-odd excited states can be written as 



N*) = Aqe Qbc / tiVitf r 2 cfr 3 exp 
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N 2 e abc 



d r\d r 2 d r 3 exp 
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(109) 



In the above formulas X x , X y and A 2 are the x, y and z components of A, respectively. 
Ni and N 2 are normalization factors of the states with Ni = 2 1 / 2 a 4 /(3 9//4 7r 3//2 ) , N 2 = 

a V(2 1 /2 3 9/4 7r 3/2 ) 
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TABLE XIII: nEDM contributed from first excited CP-odd states in the non-relativistic quark 
model, where C 4 is the Wilson coefficients of the quark models, a is defined below Eq. (|83p . The 
unit of nEDM used here is e-GeV -1 , which is different from the traditional one e-cm due to that 
the Wilson coefficients of the four-quark operators are unknown which are always in the unit of 
GeV -2 . The translation between the two units is e • GeV -1 ~ 2 x 10~ 
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e-cm. 
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TABLE XIV: Same as Table XIII, except the renormalization scale is now at ^ttF^. 



The results are shown in Table 1X111} which agree with the results extracted from Ref. 



50J. We also need to take into account the evolution of the operators between AirF n and the 
energy scale of the quark model. The results are shown in Table IXTVl with a = 0.41 GeV. 
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B. Meson Condensate Contribution through Photo-Pion Production 



In photon-pion production, there are CP-even electric-dipole couplings between the 
baryon-octet and electromagnetic fields through using f±" Some of these couplings can 
generate the neutron EDM if they violate the chiral symmetry through the quark masses and 
at the same time the meson fields acquire vacuum condensates through the CP-odd four- 
quark operators. In more physical language, the contact terms for the pion-photoproduction 
processes give rise to the neutron EDM through the diagram in Fig. HI Although the elec- 




FIG. 4: Pion-photoproduction diagram with the pion field annihilated by the four-quark operator 
into the vacuum, where the cross is a four-quark operator. 



tromagnetic field also violates chiral symmetry, it cannot generate an EDM through meson 
condensates by itself — a quark mass factor is essential. 
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TABLE XV: nEDM induced by meson condensates through pion-photoproduction. Contribution 
from operators constructed by tensor operators are neglected due to the large- Nq suppression. 



The terms of interest are made of linear products of baryon fields B and B, X- an d 

C% = ^p- 2 [«ili[B V ft{X-.ffl + 5 2 Tr[Ba^ 5 f^B]Tr[x-} + ...] , (110) 

where we have shown two of the ten possible terms. It is difficult, however, to extract the 
Wilson coefficients Si directly from experimental data. Some of the coefficients have been 
estimated by calculating the contribution from the excited baryon states in the context 
of the two-flavor scenario [53]. In the two-flavor scenario, neglecting the isospin violation 
generated by the difference between the up and down quark masses, the terms relevant to 
nEDM can be written as 

^-flavor = [(< _ ^ f ,u + ^ (/f )])( _ N ? (m) 
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where N = ( ), and in the two-flavor case, f+ u = e(^Q£ + £,Q£})F fJ,l/ , in which Q = 
(1 + t 3 )/2. Expanding and X-i we can get the nEDM induced by the condensate of 7r°; 

d = Sea^B (m u + m d )(n°) 



From Ref. |53J, one can get the contribution to a± from A and p internal states, which is 

ai = -0.156GeV~ 3 . (113) 
Using this, one can estimate the nEDM induced by the pion condensate, as shown in Table 

KYI 

C. CP-Odd Baryon Mass Contribution 

The CP-odd baryon-mass terms considered in the previous section generate a CP-odd 
part of the baryon wave function. This part can transform a magnetic moment term into 
an EDM contribution. The physics of this is shown in Fig. [5j 



(a) (b) 

FIG. 5: The CP-odd mass of neutron turns the tree level magnetic moment into an EDM. The 
cross is the tree level magnetic moment, the gray dot is the CP-odd mass of the neutron and the 
black dot is the CP-odd pion-nucleon coupling. 



The mass terms of the neutron can be written as 

£mass = -m n nn - m' n nij 5 n . (114) 

Note that the neutron field n here is already redefined using the transformation in Eq. (|36|) 
after taking into account the meson condensate effect as discussed in the previous section. 
Redefining the neutron field again through a chiral rotation, 

n = exp(-i^-)7 5 n , (115) 

the mass term becomes the standard one, 

£mass = -m n n'n' . (116) 

On the other hand, the tree level anomalous magnetic moment of the neutron can be written 

as 

1 k 



r = -ncr^nF (117) 

mag. mom. . Ibu lb± ^iv ■ v-"-- 1 - 1 / 
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The redefinition in Eq. (11151) generates a neutron EDM, 



jEDM _ K n m n M-\o\ 

"CP-oddmass ~ 2m 2 ' V LL °) 

The experimental values of the anomalous magnetic dipole moments of the nucleons are 
k p = 1.7928, K n = —1.9131. The numerical values of this contribution have been shown in 
Tables IXVII and IX VIII For the tensor operator ui^sa^uda^d, the contribution from the 
CP-odd mass of the nucleon is particularly large. The CP-odd mass also gets a quantum 
correction shown in diagram (b) of Fig. |5j It is easy to see that this term does not have any 
chiral enhancement and is of a higher-order effect. 



D. Leading Chiral Loop Contribution 

The contribution we have considered so far has a smooth chiral limit, i.e., regular as the 
quark masses go to zero. The leading contribution in the chiral limit, however, involves the 
pion loop with an infrared divergence. This contribution was first calculated by Crewther 
et al [15J, and has been studied thoroughly in the literature (see Fig. E]). Diagrams (a) and 
(b) in Fig. E] contain an infrared divergence which is regularized by the mass of pion and an 
analytical part. The constant part is canceled by diagrams (c) and (d). Diagrams (e) and 



(f) cancel with each other [13J. Therefore, up to terms of order (m w /m n ), the neutron EDM 



n p <^p n n P <^P » 

(c) (d) 



P n n <3 P 

(e) a 



FIG. 6: Charged-pion loop contribution to neutron EDM (without the anomalous magnetic mo- 
ment), where the black dots represent the CP-odd vertices. 



generated by the charged pion loop can be written as 13 

eV2 



dt, 



16n 2 F 7 



■h c {D + F) In {ml/m 2 n ) 



(119) 



where D + F = —gA = —1.26 is the CP-even pion-nucleon coupling (the signs of D and F 



is different from that in Ref. [24| because we are using a different definition of the chiral 
transformation of U), and h c is the CP-odd pion-nucleon coupling defined in Eq. ( I7H|) . Note 



that, in Fig. [61 the contribution from the proton's anomalous magnetic moment has not 
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FIG. 7: Contribution from the tree level anomalous magnetic moments of proton and neutron, 
where the crosses are anomalous magnetic moments of nucleons and the dots are CP-odd vertices. 



been included. To include this contribution, we consider all these diagrams in Fig. [7] where 



the neutral pion loop is also present, and the result is Ref. [21 

e D + F 



d. 



^ 167T 2 F 7 



- — y-V2h c K p + h n K n j F n ( — | J 



(120) 



where n n and K n are tree-level anomalous magnetic moments of protons and neutrons, re- 
spectively, and 



, * 3 3s - s 2 , s(5s - s 2 ) - 4s Js-s 2 /A 

r n (s) = s In s H — arctan — ; . 

ny ' 2 2 2^s - s 2 /4 s/2 

We can see that there is no chiral enhancement in Fr,(s). 



;i2i) 
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TABLE XVI: nEDM from the P-odd and CP-odd four-quark operators composed of color-singlet 
currents. Different contributions are shown. 
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TABLE XVII: Neutron EDM generated by P-odd and CP-odd four-quark operators composed of 
color-octet currents. The labels have the same meaning as in Table IXVll 



Using the above, we estimate the pion-loop and the CP-odd mass contributions to neutron 
EDM due to the P-odd and CP-odd four-quark operators. The results are listed in Tables 
IXVII and IXVIII Although the charged pion-loop (Fig. [6]) dominates in the chiral limit, its 
numerical value is actually about an order of magnitude smaller than the analytical chiral- 
loop contribution (Fig. [7]). This is due to the enhancement of h n relative to h c in the large 
Af c limit. 

The P-odd and CP-odd four-quark operators can also lead to nonvanishing P-odd and 
CP-odd interaction like n — > KT, and n — > Arj. These interactions can generate nEDM 
through kaon- or eta-loop diagrams. However, there is no reason to believe that the kaon- or 
eta-loop contribution should be more important than the pion-loop contribution so that it 
would not change the order-of-magnitude estimate of nEDM generated by those four-quark 
operators without the strange quark. For those operators containing strange quark the 
estimation may not be reliable and the kaon- or eta-loop contributions should be included. 



E. Comparison with Other Calculations and the Error-bars of this Calculation 

The P-odd and CP-odd four-quark contributions to neutron EDM have been studied using 



different approximation methods in the literature [19H21I]. The problem is that it is difficult 
to get an estimate on the errors in any of these methods. This is the strong motivation for 
the alternative study presented here. By using a completely different approach, we hope to 
get a better idea how well one actually estimates these hadronic matrix elements. 

In Ref. jlji], the authors used the external field method, factorization and QCD sum 
rules to make a direct calculation of the neutron EDM. Their result is supposed to be the 
total contribution, although it is unclear how the chiral physics would be included in this 
approach. Their numbers are listed in Table IXVIIII as "factorization and QCD sum rule." 
The result is, in general, comparable to the charged pion-loop contribution, although the 
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contribution to the tensor operator is particularly large. 

In Ref. [20[, the authors also calculated the contributions of the pion-loop as we do in this 
paper. They used entirely the factorization method to calculate the CP-odd pion-nucleon 
couplings, including the effects that the CP-odd operators can annihilate the neutral pion 
in the vacuum. Taking the operator ui'j^udd as an example, their factorization works like 
this: 

(nn°\ui'~f5udd\n) = (n\dd\n)(ii \ui'y5u\0) 

+ (0|dd|0) (imr°\ui^ 5 u\n) - ^{n^\C^ CD \n^){^\ui^u\{M ,(122) 

where jCqcd is the usual QCD Lagrangian. The terms inside the bracket on the second line 
of the above formula cancel each other. The reason is that uij$u is just a CP-odd mass of the 
up-quark which can be rotated away through chiral transformation, except for a possible 
Ua(1) contribution. Thus these two contributions should cancel with each other exactly. 



This is first noticed in Ref. [54J in the spirit of the Feinberg-Weinberg-Kabir theorem |55| . 
Using this method, one can get the CP-odd vertices, h c and h n as shown in Table IVII For 
the charged coupling h c , one needs to do a Fierz transformation, from which one can get a 
suppression factor of 1/12, where 1/3 is from the color factor and the other 1/4 is from the 
spin. Therefore, h c is one order of magnitude smaller than h n . The corresponding nEDM 
calculated using this method is included in Table IXVlHI as well. 

From Table IXVIIH taking the operator ui^^udd as an example, one can see that the 
magnitude of our result is comparable with what obtained using naive factorization method 
but with a different sign; also the our result is about one order of magnitude larger than 
the result estimated using QCD sum rules. In our calculation, we separate the contribution 
into the meson condensate contribution and the direct matching contribution. The vacuum 
saturation method is used to calculate the meson condensate contribution to h c and h n . 
This vacuum saturation method using to calculate the meson matrix elements is accurate 
in the large- A*^ limit, which means the calculation for this contribution is accurate up to 
1/Nc (39[. From Table IXVIj one can see that the meson condensate contributions dominate 
over the direct matching contributions. Therefore, for operators generating unsuppressed 
meson condensates (see Sec IV for detailed discussions), a conservative uncertainty can be 
set to be a factor of two. 



In Ref. [20(, the authors also used the vacuum saturation approach to get the factorization 
result as shown in Eq. f ll22p . However, in the case of baryon matrix element, the non- 
factorized contribution is not suppressed in the large- Nc limit [39], therefore the missed 
non-factorized contribution should be of the same order as the factorized contribution shown 



in Eq. (I122p . The calculation using QCD sum rules in Ref. [19| did not include the meson 
condensate contribution, therefore, their calculation might miss an important contribution. 
The factor of two uncertainty can also be seen from the Feinberg-Weinberg-Kabir the- 



orem 



55| . Applying to this context, the theorem dictates that CP-odd (3,3) two-quark 
operators give no contribution to CP-odd processes. However, since we are using a hy- 
brid method, this theorem may not be satisfied. Therefore, the amount of violation of this 
theorem can be seen as an estimate of the error of this calculation. Take the operator 
uijsu — dij 5 d as an example, following the prescription in Sees. Ill and IV, one can get 
meson-condensate contribution to the neutral CP-odd pion-nucleon coupling which can be 
written as 

Ksr = ^±^„_i^ i (123) 
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TABLE XVIII: Comparison of different methods, nEDM calculated by factorization in Ref. 2(], 21] 
are shown as "naive factorization". The column on the right side shows nEDM calculated using 
factorization and QCD sum rules [19]. The unit of the numbers is 10 _3 eC4GeV. 



where m„ = = m is assumed for the sake of simplicity, C3 is the Wilson coefficient 
of the two-quark operator and the definitions of c\ and C3 can be found in Eq. ( I76i) . If 
the cx-term is also employed to do the direct matching, one can easily show that the direct 
matching contribution cancels the meson condensate contribution exactly. Instead, in order 
to get the uncertainty of our calculation we need to do the direct matching using the quark 
model. Since the operator includes only products of two quark fields, the calculation using 
the quark model is straightforward, which gives 

C = - ^ , (124) 

£ 7T TV 

where the factor of 3 is due to that in the quark model the nucleon contains three con- 
stituent quarks. R ~ 1.7 comes from the perturbative QCD effect as discussed in Sec. V. 
The anomalous dimensions of the operator discussing here is as the same as the anomalous 
dimensions of the quark mass. The relative sign between the direct contribution and the 
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meson condensate contribution is as desired. However, the magnitude of the direct contri- 
bution is about two times smaller than the meson condensate contribution. The mismatch 
between the two contributions is due to that quark model does not differentiate (N\qq\N) 
and (N\q^q\N). From this mismatch one can see that the inaccuracy of the direct contribu- 
tion calculated using quark model might be a factor of two. Therefore, conservatively, the 
total inaccuracy for those operators having unsuppressed vacuum condensate contributions 
can be seen as a factor of two. 
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TABLE XIX: Upper bound on the Wilson coefficients of P-odd, CP-odd four-quark operators, 
calculated using the experimental data and hadronic matrix elements in this work. 



VI. CONCLUSION 

In this paper, we studied the four-quark contributions to the neutron EDM, which dom- 
inate over other QCD operators in some new physics models. Our approach was based on 
chiral expansion and simple quark models. It is well known in the literature that the leading 
chiral contribution comes from one-pion loop which dominates in the chiral limit m n — > 0, 
just like in the case of the nucleon electric polarizability. Therefore, one needs to calculate 
the four-quark contribution to the CP-odd pion nucleon couplings. We studied these cou- 
plings in simple quark models, as an alternate to large- N c factorization. We also considered 
0(1) contribution from direct matching and pion-condensation to the dipole moment, as 
well as the CP-odd nucleon mass contribution through the magnetic moment. The resulting 
11EDM can be compared with those from the naive factorization and QCD sum rules. The 
comparison provides us some idea on the hadronic physics uncertainty in the neutron EDM 
calculation. Our approach also provides a formalism for lattice QCD calculations of the 
nucleon matrix elements of the four-quark operators. 

Using the matrix elements thus obtained, we obtain new-physics-independent upper 
bounds on the Wilson coefficients of four-quark operators from the experimental data. The 
current experimental upper bound on neutron EDM is 2.9 x 10~ 26 e cm j30[. If we assume 
that there is no significant cancelations among the contributions from these operators, we 
can use the experiment limit to give upper bounds to the Wilson coefficients of individual 



38 



operators. In our calculation, the strange quark effects were ignored, and we considered only 
operators composed of up and down quarks. The final results are shown in Table IXIXl 

It is interesting to note that the chiral-enhanced contribution is actually large-Y c sup- 
pressed. In fact, the non-singular part of the chiral-loop contribution numerically dominates 
over the singular one. This suggests a large- N c analysis of the neutron EDM, including the 
delta resonance contribution. However, this is beyond the scope of this paper. 
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